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Abstract 

Cosmic-ray proton and antiproton spectra were measured at mountain altitude, 
2770 m above sea level. We observed more than 2 x 10 5 protons and 10 2 antiprotons 
in a kinetic energy range between 0.25 and 3.3 GeV. The zenith-angle dependence 
of proton flux was obtained. The observed spectra were compared with theoretical 
predictions. 
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1 Introduction 

Primary cosmic rays hit the Earth's atmosphere and produce baryons and 
mesons via hadronic interactions. Absolute fluxes of these secondary cosmic 
rays at any altitudes can be calculated based on primary cosmic-ray inten- 
sity and interaction cross sections. Observations of the secondary cosmic rays 
are very important for understanding the production and propagation of sec- 
ondary cosmic-ray particles inside the atmosphere. It will provide information 
to verify or to improve the theoretical models. A number of experiments for 
measuring primary cosmic rays have been carried out at balloon altitude or 
in space, and atmospheric muons on the ground as well as under the ground. 
A long-period observation at mountain altitude with high statistics will pro- 
vide a good reference for a study of secondary cosmic-ray particles inside the 
atmosphere. 

Only a few measurements of cosmic-ray proton fluxes at mountain altitude 
were reported [1,2,3,4]. Antiproton flux at mountain altitude was also reported 
[4] previous to this work. The reported flux was quite higher than theoretical 
predictions. 

We report here a new observation of cosmic-ray protons and antiprotons at 
mountain altitude in a kinetic energy range between 0.25 and 3.3 GeV. A 
spectral shape of cosmic-ray antiprotons at mountain altitude was measured 
for the first time. These data will provide essentially important information 
about hadronic interactions between cosmic rays and atomic nuclei inside the 
atmosphere. 



2 BESS Experiment 

2.1 Detector 

The BESS (Balloon-borne Experiment with a Superconducting Spectrometer) 
detector [5,6,7,8,9,10] is a high-resolution spectrometer with a large acceptance 
to perform highly sensitive searches for rare cosmic-ray components, as well 
as precise measurement of the absolute fluxes of various cosmic rays. Fig. 1 
shows a schematic cross-sectional view of the BESS instrument. In the cen- 
tral region, a uniform magnetic field of 1 Tesla is provided by using a thin 
superconducting solenoidal coil. A magnetic-rigidity (R = Pc/Ze) of an in- 
coming charged particle is measured by a tracking system, which consists of a 
jet-type drift chamber (JET) and two inner-drift-chambers (IDC's) inside the 
magnetic field. The deflection (-R -1 ) and its error are calculated for each event 
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by applying a circular fitting using up to 28 hit points, each with a spatial 
resolution of 200 /zm. The maximum detectable rigidity (MDR) was estimated 
to be 200 GV. Time-of-flight (TOF) hodoscopes provide the velocity (f3) and 
energy loss (dE/ dx) measurement. A 1/(3 resolution of 1.6 % was achieved 
in this experiment. In order to separate protons and antiprotons from muons, 
the BESS spectrometer is equipped with a threshold-type aerogel Cerenkov 
counter. The refractive index of silica aerogel radiator is 1.022, which corre- 
sponds to a threshold kinetic energy of 3.6 GeV for protons and antiprotons. 
Each particle is identified by requiring proper 1/(3, as well as dE/ dx, as a 
function of the rigidity. 

There are a plastic scintillating counter, an acrylic Cerenkov counter and a 
lead plate just above the bottom TOF hodoscope. The lead plate covers about 
1/5 of the total acceptance. These counters were used for electron/muon iden- 
tification in the muon analysis [11]. In the analysis of protons and antiprotons 
discussed here, the signals from these counters were not examined. 

All detector components are arranged in a cylindrical configuration. This sim- 
ple configuration and the uniform magnetic field result in a large geometrical 
acceptance and uniform performance in momentum measurement. 

A simple coincidence of the top- and bottom-TOF hodoscope signals issues 
the first-level "TO-trigger." The live data-taking time is measured exactly by 
counting 1 MHz clock pulses with a scaler system gated by a "ready" status 
that controls the TO-trigger. 

2.2 Observations 

The atmospheric cosmic-ray events were observed at Norikura Observatory, 
ICRR, the University of Tokyo, Japan. It is located at 36° 06' N, 137° 33' E. The 
altitude is 2,770 m above sea level. The vertical geomagnetic cutoff rigidity is 
11.2 GV [12]. 

The observation was performed during two periods of 17th - 19th and 21st 
- 23rd of September 1999. During the observation, the atmospheric depth 
and temperature varied as shown in Fig. 2. The mean (root-mean-square) 
atmospheric depth and temperature were 742.4 (2.9) g/cm 2 and 10.9 (1.1) °C, 
respectively. The BESS detector was operated outside the building and in a 
nylon plastic sheet tent so as to be less influenced by ambient conditions. 
There was no thick material within a field of view of the detector. A gas flow 
system was implemented to keep the purity and pressure of the gas inside the 
chambers stable enough for operation. 

The TO-trigger rate was about 50 Hz. It is much lower than that at a usual 
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balloon floating altitude of about 37 km. Thus, all events which satisfied the 
TO-trigger condition were recorded without imposing second-level trigger nor 
on-line rejections. This simple trigger condition reduced systematic errors. The 
live-time ratio was as high as 98.8 % through the observation. The number of 
collected events was about 2.0 x 10 7 . 



3 Data analysis 

Particle mass was reconstructed by the velocity, rigidity and charge for each 
event, and was required to be consistent with a proton and an antiproton. 
Fig. 3 shows 1//3 distribution as a function of rigidity after requiring that there 
was no light output from the aerogel Cerenkov counter, so as to discriminate 
protons and antiprotons from muons. 

The procedure of data analysis was almost the same as that of the balloon- 
flight data [13,14,15] except for a method of evaluating the performance of 
the aerogel Cerenkov counter. The cosmic rays at the balloon altitude are 
dominated by protons. On the other hand, those at mountain altitude con- 
sist almost of muons. Therefore the efficiency of the aerogel Cerenkov counter 
for protons and antiprotons was not obtained from the observed data in this 
analysis. The efficiency and its dependence on particle's velocity were esti- 
mated by using muons in a lower rigidity range of \R\ < 1.7 GV, where muons 
were identified without contaminations. The obtained efficiency was applied to 
protons and antiprotons. This process arose an additional systematic error of 

4 %. In a higher rigidity range of \R\ > 1.7 GV, protons and antiprotons were 
contaminated by muons. The number of contaminating muons was estimated 
by multiplying the total number of observed particles by the probability that 
a muon fakes a proton or an antiproton. The fake probability was evaluated to 
be 10~ 4 by using energetic particles above 10 GeV and the same muon sample 
as was used for estimating the efficiency of the aerogel Cerenkov counter. The 
muon background amounted to about 20 % for the antiprotons above 2 GeV. 
It was negligibly small for the protons. 

In this analysis, the zenith angle (8 Z ) was limited as cos0 2 > 0.95 for protons. 
Because of the very small flux of antiprotons, the lower limit of cos6 z for 
antiproton analysis was relaxed to be 0.84 so as to improve its statistics. 



4 Results and discussion 

The observed energy spectra of protons and antiprotons are shown in Figs. 4 
and 5, respectively, together with previous measurements at mountain altitude 
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[1,2,3,4] and theoretical predictions [16,17,18,19]. Tables 1 and 2 summarize 
the obtained spectra. The first and second errors in Tables 1 and 2 represent 
statistical and systematic errors, respectively. Since the statistics of antiproton 
measurement was very limited, the statistical error was calculated as a 68.7 % 
confidence interval based on a "unified approach" [20]. 

4-1 Protons 

There is some difference among the proton fluxes observed by various exper- 
iments as shown in Fig. 4. Table 3 summarizes the locations of cosmic-ray 
observations at mountain altitude. According to our simple Monte Carlo sim- 
ulations, the difference was generally explained by the different altitudes and 
cutoff rigidities. The overall spectral shape is reproduced by theoretical calcu- 
lations [16,19]. 

Fig. 6 shows the zenith angle dependence of the observed proton flux, F(cos 9 Z ), 
in two kinetic energy regions. In a simple one-dimensional approximation, the 
zenith angle dependence can be expected as 



where X is the atmospheric depth at the observation site and A is the atten- 
uation length of protons inside the atmosphere. The dashed lines in Fig. 6 
show the expectation, in which X/A = 6 was assumed, since the atmo- 
spheric depth was 742.4 g/cm 2 and the proton attenuation length is about 
120 g/cm 2 . Normalization was done by tuning F so that the total flux, 
2n j F(cos9) cos 9d cos 6, is equal to the observed total flux. Relatively good 
agreement is found between the observed data and the expectation in the 
higher energy region. In the lower energy region, however, much larger devi- 
ation is observed in the zenith angle dependence. These facts are most likely 
due to the effect of angular spread of secondary protons produced via nu- 
clear interactions. This "three-dimensional effect" was evaluated by a simple 
analytic method considering the distribution of transverse momentum at the 
production point of secondary protons. The solid lines shown in Fig. 6 give 
the results of the analytic calculation, in which the angular spread was taken 
into account. They reproduce the observed data in the whole energy range 
better than the expectation by the one-dimensional calculation. 

As mentioned above, the proton flux shown in Fig. 4 was obtained by using 
the events satisfying the zenith angle condition of cos #2 > 0.95. The average 
zenith angle, (cos 6^), was 0.98. The vertical flux F = F(cos9 z — > 1) was 
obtained by fitting the observed fluxes using a formula (1), where "X/A" was 




(1) 
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used as a fitting parameter. The F obtained by fitting the data in cos 9 Z > 0.95 
is summarized in the last column of Table 1. The systematic error due to this 
fitting procedure was estimated to be around 5 % by checking variations of 
F in changing the fitting regions. 



4-2 Antiprotons 

Fig. 5 shows that the obtained antiproton spectrum above 1 GeV is generally 
consistent with theoretical predictions obtained through transport equation 
calculations by Bowen and Moats [16] and Stephens [17], as well as through a 
three-dimensional Monte Carlo simulation by Huang et al. [18]. In the lower 
energy region, however, the transport equation calculations show significant 
disagreement with our results. Those three calculations predict similar an- 
tiproton spectrum above 0.3 GeV at a thin residual atmosphere of 5 g/cm 2 , 
but their predictions do not agree so well with each other at mountain al- 
titude [16,17,18]. It means that the shapes of antiproton production spectra 
are very similar to each other, and that the difference at mountain altitude is 
to be made mainly due to the different treatment of the propagation inside 
the atmosphere. The production spectrum of antiprotons has a sharp peak 
around 2 GeV [17]. The flux at 0.5 GeV is as small as 1/10 of the peak flux. 
Therefore most antiprotons observed below 1 GeV are tertiary antiprotons, 
those which have been produced inside the atmosphere and then lost their 
energy during the propagation inside the atmosphere. In this case, interac- 
tion processes of p + A (nuclei) have to be precisely treated for an accurate 
evaluation of antiproton spectrum at mountain altitude. In the calculation 
of Refs. [16,17], the probability that an antiproton with initial energy of E 
possesses energy E after a collision was assumed to be uniform from E = to 
E = 0.90E ° r 0.95i?o- On this assumption, the average energy after a collision 
is about a half of the initial energy. In the calculation of Ref. [18], on the other 
hand, only annihilation channels were taken into account in the inelastic in- 
teractions. Non-annihilating inelastic processes were not included as a process 
of antiproton energy loss, thus a tertiary antiproton remains as energetic as 
before the interaction. Our result suggests that the non-annihilation process 
does not make significant contribution to the inelastic interactions of antipro- 
tons during their propagation inside the atmosphere. Thus, the antiprotons 
are considered not to lose their energy considerably in the non-annihilating 
inelastic processes. 

The antiproton flux measured by Sembroski et al. [4] is much higher than our 
measurements in spite of the similar observational conditions such as altitude 
(Table 3). In that experiment, a particle trajectory was measured with six 
chambers only outside the magnetic field. No information about the parti- 
cle trajectory was obtained inside the magnetic field. It might lead to some 
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difficulty in strict identification of antiprotons. 



5 Summary 

We have measured proton and antiproton spectra in a kinetic energy range 
of 0.25 - 3.5 GeV at Mt. Norikura, 2,770 m above sea level, in Japan. The 
vertical geomagnetic cutoff rigidity is 11.2 GV. The mean atmospheric depth 
during the measurement was 742.4 g/cm 2 . 

The measured proton spectrum is consistent with the previous results, if the 
different altitudes and cutoff rigidities at their observation sites are taken 
into account. The zenith angle dependence of the proton flux was observed. 
It suggests importance of a three-dimensional effect of the angular spread in 
secondary baryon productions especially in a low energy region. 

The measured antiproton flux above 1 GeV generally agrees with the theo- 
retical predictions. In the lower energy region, however, our flux shows much 
better agreement with the calculation performed in assuming that tertiary an- 
tiprotons remain as energetic as before interactions. Our result suggests that 
the energy loss of antiprotons due to non-annihilation process is not significant 
during their propagation inside the atmosphere. 
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Fig. 1. Schematic cross-sectional view of the BESS detector. 
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Fig. 2. Atmospheric depth and temperature during the observation. The experiment 
was performed during two periods of 17th - 19th and 21st - 23rd. 
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Fig. 3. Identification of antiproton events. The solid lines define the antiproton 
selection band. 
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Fig. 4. The observed proton spectrum compared with the previous works, and with 
the theoretical calculations. 
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Fig. 5. The observed antiproton spectrum compared with the previous works, and 
with the theoretical calculations. The dashed line was reproduced by interpolating 
between two spectra at 700 g/cm 2 and 900 g/cm 2 using a power law. 



13 




Fig. 6. The zenith angle dependence of proton flux. The dotted and solid lines 
show the expected dependence in simple one-dimensional and three-dimensional 
calculations, respectively. 
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Table 1 

Observed spectrum of protons. 
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Table 2 

Observed spectrum of antiprotons. 
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Table 3 



Summary of the observation sites at mountain altitude. 
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Refs. [1,2] 


Refs. [3,4] 


Site 


Mt. Norikura 


Mt. Aragats 


Mt. Lemmon 


Altitude 


2770 m 


3200 m 


2750 m 


Atmospheric Depth 


742 g/cm 2 


710 g/cm 2 


747 g/cm 2 


Cutoff Rigidity 


11.2 GV 


7.6 GV 


5.6 GV 
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